Universal Quantum Entanglement between an Oscillator and Continuous Fields 
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Quantum entanglement has been actively sought for in optomechanical and electromechanical 
systems. The simplest such system is a mechanical oscillator interacting with a coherent beam, 
while the oscillator also suffers from thermal decoherence. For this system, we show that quantum 
entanglement is always present between the oscillator and continuous outgoing fields — even when 
the environmental temperature is high and the oscillator is highly classical. Such universal entangle- 
ment is also shown to be able to survive more than one oscillation cycle if characteristic frequency 
of the optomechanical interaction is larger than that of the thermal noise. Furthermore, we derive 
the effective optical mode that is maximally entangled with the oscillator, which will be useful for 
future quantum computing and encoding information into mechanical degrees of freedom. 

PACS numbers: 03.67.Bg, 03.65.Ta, 42.50.Wk 
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Introduction. Entanglement, as one of the most fascinat- 
ing features of quantum mechanics, lies in the heart of 
quantum computing and many quantum communication 
protocols [l[ . Great efforts have been devoted to theoret- 
ical and experimental investigations of quantum entan- 
glements in different systems with discrete or continuous 
variables. Due to advancements in fabricating low- loss 
optical elements and high-Q mechanical resonators, the 
quantum entanglement in optomechanical systems has 
recently aroused great interests, especially when many 
table-top experiments demonstrated significant cooling of 
mechanical degrees of freedom via active feedback or pas- 
sive damping (self-cooling) @, i, i, @, 0, 0, 0, E3] , which 
unveils the possibility of achieving quantum ground state 
of macroscopic objects llj, [TJ, ll3|- This not only paves 
the way for high-precision measurements but also in- 
corporating mechanical degrees of freedom as possible 
medium for storing and retrieving quantum information. 

Theoretical analysis shows that by coupling oscilla- 
tors to a Fabry-Pcrot cavity, one can create stationary 
(Einstein-Podosky-Rosen) EPR-type quantum entangle- 
ment between optical modes and an oscillator 1J] or even 



la lfl. InRef. m 



between two macroscopic oscillators 
it was shown that entanglement between two oscillators 
can also be created by conditioning on continuous mea- 
surements of the common and differential optical modes 
in a laser interferometer. Interestingly, such entangle- 
ment does not depend on the environmental tempera- 
ture T explicitly but rather scales as the ratio between 
characteristic interaction frequency Q q (equivalent to op- 
tical power), and characteristic thermal- noise frequency 
Clp (equivalent to T). In contrast, T enters explicitly 
and entanglement generally vanishes at high temperature 
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in cases considered in Refs 
ancy arises from the following facts 
Ref. 17j is tuned and thus always stable 



This discrep- 
(i) The cavity in 
while in Refs. 




FIG. 1: A schematic plot of the model and the corresponding 
spacetime diagram. Here x, a, and b denote the oscillator po- 
sition, ingoing and outgoing fields respectively. For clarity, we 
intentionally place a and b on difference sides of the oscillator 
world line. The tilted lines represent the light rays. Up to 
some instant we are concerned with (t = 0), the optical fields 
entering later are out of causal contact and thus irrelevant. 
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16j . stability requirements of the system set an 
upper limit on fl q ] (ii) More importantly, due to finite 
transmission of the cavity, information leaks into the en- 
vironment. Therefore, even regardless of thermal heat 
bath, the reduced system consisting of cavity modes and 
the oscillator is not in a pure state. In Ref. 17] , however, 
there can be either no cavity or cavity with very broad 
bandwidth, the outgoing fields containing information of 
oscillator motion are all registered by photodetector. 
This motivates us to consider that the entanglement 



scaling obtained in Ref. 17j can be inherent in the sim- 
plest system — an oscillator interacting with a coherent 
beam, which models the essential process in all above- 
mentioned optomechanical systems. The model and its 
spacetime diagram are shown schematically in Fig. [T] 
Similar system was analyzed previously by Pirandola et 
al. [3|. They used narrow-detection-band approxima- 
tion to introduce sideband modes, which maps outgoing 
fields into two effective degrees of freedom. In the situa- 
tion here, sideband modes are not well-defined, because 
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the interaction turns off at t — and only half-space 
[—00, 0] is involved. Instead, we will directly evaluate the 
entanglement between the oscillator and outgoing fields b 
(infinite degrees of freedom) using the pos itivity of partial 
transpose (PPT) criterion fljj Hfl ED, M 0, [Ij \m [H] . 
Only in weak-interaction and low-thcrmal-noisc limit 
(fl q , rip <C w m ) can wc make correspondences between 
our results and those obtained in Ref. [l8j |. 
Dynamics and Covariance Matrix. The Heisenberg equa- 
tions for this optomechanical system are simply 

x(t)=p(t)/m, (1) 
p(t) = -2"f m p(t) -mu 2 m x(t) +aa l (t) +f t h(*), (2) 
61 (t) = ai (t) , 6 2 (f ) = a 2 (t) + (a/K) x(t) . (3) 

Here x and p are oscillator position and momentum; ai 
and 6,; (i = 1, 2) are quadratures of ingoing and outgoing 
optical fields, a\ = (a + aJ)/y/2 and S2 = (a — a))/ (i\/2) 
(the same for 61,2); a = (IqHujq/c 2 ) 1 ^ 2 = (hmVt 2 ) 1 / 2 is 
the optomechanical coupling strength, where loq and Iq 
denote the laser frequency and optical power respectively 
and we have defined f2 g ; a Si is the radiation-pressure 
term. Since [ai(t), ai(t')] = 0, the presence of thermal 
noise £th ensures the correct commutator between x(t) 
and p(t) [27| . The solution to oscillator position x is 

*(*) = J t _ oo dl/G x (t-t')[aa 1 (t')+i th (t')}, (4) 



where Green's function G x (t) 



= P ~Jmt 



sm.(uj m t)/(muj m ). 



The radiation-pressure term a a\ induces quantum corre- 
lations between the oscillator and the optical fields, but 
it is undermined by £ t h- The question would be whether 
quantum entanglement exists or not after evolving the 
entire system from t = —00 to 0. Since variables involved 
are Gaussian and linear dynamics will preserve Gaussian- 
ity, entanglement is completely encoded in the covariance 
matrix V of the optomechanical system. With optical 
fields labeled by continuous coordinate t, elements of V 
involving optical degrees of freedom would be defined in 



the functional space £ 



V 



X), 0]. Specifically, 

A C T 
C B 



(5) 



i.e. gjijy > 0. In the phase space of continuous Gaussian 
variables, this reduces to the Uncertainty Principle 



V pt + (1/2)K > 0. 



(6) 



Here commutator matrix K = 0)^=1 2 &y with o~ y denot- 
ing Pauli matrix. According to the Williamson theorem, 
there exists a symplectic transformation S G S , p (2jv+2,r) 
such that S T V pt S = @^=\ Diag[Afc, A&]. Using the fact 
that S T KS = K, the above Uncertainty Principle reads 
Ajt > 1. If this fails to be the case, i.e. 3A& < 1, the 
states are entangled. The amount of entanglement can 
be quantified by the logarithmic negativity Ej\f [2~1 ] and 



£V = max[— ^ fe lnAfe, 0] for A: : Afe < 1. 



(7) 



In the case here, N approaches 00. Besides, the partial 
transpose is equivalent to time reversal. Therefore V pt = 
V|p(o)_ >-p(o)- Normalizing x and p with respect to their 
zero-point values, the commutator reads [x,p] ~2i. For 
the optical fields, we set [b\(t), &2(0] = 2i 5(t — t'), which 
gives the coordinate representation of K. 

According to Ref. [24j , A& can be obtained by solving 

V pt v = (l/2)AKv, (8) 

where v = [a>o, /3o, \a), \(3)] T with |/) denoting the vec- 
tor in C 2 [— 00, 0]. Due to uniqueness of |a) and |/3 X 
terms of ao and (3q for any A < 1 (non-singular), Eq. 
leads to the following characteristic equation 



m 



det[A + X(T y — C T (A a y + H)~ 1 C] =0 



It can be shown that 



(AfTy + B) _1 = 



b{m~ 



-b{m- 

M- 1 



(9) 



(10) 



where we have used the fact that B\ 2 = B21 in C 2 [— 00, 0] 
and have defined B\ = Bn — i A and M = B22 — Bt B\. 
The integral operator M can be inverted via Wiener- 
Hopf method. Given any function \g) = M _1 |/i), in the 
frequency domain, it reads 



J^dte^M-^h) = [l/j,-][h/j>+]-. (11) 



Here A* 



(Xi Xj) sym (i,j 



1,2) with vector X = 
[x(0), p(0)] and (Xi A,) sym = (X t X 3 + X 3 X t )/2 denot- 
ing symmetrized ensemble average; Cjj and By should 
be viewed as vectors and operators in £ 2 [— 00, 0]. In the 
coordinate representation, (t\ Gy) = (Xibj(t)) sym and 
(t| Bj 3 \t') = (bi(t)bj(t')) B y m , in which (|) denotes the 
scalar inner product in £ 2 [— 00, 0]. 

PPT Criterion. According to Refs. 23|, 26[, in order for 
one particle and a joint system of arbitrarily large N par- 
ticles to be separable, a necessary and sufficient condition 
is that partially transposed density matrix g^ N (with re- 
spect to the first particle) should be positive semidefinite, 



Here [ ] _ means taking the causal part of given function 
(with poles in lower-half complex plane) and factorization 

= A + i\(a 2 /h)(G x -G* x ) + (a/h) 2 S F G x G* (12) 

with A = 1 — A 2 and G x denoting the Fourier transfor- 
mation of G x (t). In the above equation, ip + (V 1 — ) an d 
its inverse are analytic in upper-half (lower-half) com- 
plex plane, 0) = V>2.(fi) = 'Z'-(^)- I 11 deriving Eq. 
([12]) . we have used (aj(i) aj(t')) 



sym 



Sij 5(t — t'), and for 
thermal noise, Markovian approximation is applied and 
(lth(*Kth(*')) sy ni = S F 6(t-t') with Sf = 4mj m k B T = 
2hm flp and flF defining the characteristic frequency. 
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FIG. 2: Logarithmic negativity E^f as a function of the ratio 
Qq/Qf- A mechanical quality factor Q m = 10 3 is chosen. 

Universal Entanglement. Finally, an implicit polynomial 
equation for the simplectic eigenvalue A is derived from 
Eq. ([9]). As it turns out, there always exists one eigen- 
value A that is smaller than one and it only depends on 
the ratio between f2 g and Q f , which clearly indicates the 
universality of the quantum entanglement. In Fig. [2l 
the corresponding logarithmic negativity (c.f. Eq. ([7])) 
is shown as a function of tlg/Slp- For a high-Q oscillator 
Q m = u m /(2~f m ) S> 1, up to the leading order of 1/Q m , 
a very elegant expression for Ej^ is derived and it is 



£jV = (l/2)]n[l + (25/8)n*/fi? 



(13) 



Thermal Decoherence. To investigate how long such en- 
tanglement can survive under thermal decoherence, af- 
ter turning off the optomechanical coupling at t = 0, 
the mechanical oscillator freely evolves for a finite du- 
ration t, driven only by thermal noise. Due to thermal 
decoherence, entanglement will gradually vanish. Math- 
ematically, the simplectic eigenvalue will become larger 
than unity when r is larger than the survival time t s . 
By replacing [&(0),p(0)] with [x(t),p(t)} and making 
similar analysis, up to the leading order of l/Q m , t s 
satisfies a transcendental equation: 4 Cl F 9 2 — (2 il F + 
tt 2 ) 2 sin 2 ft. 
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0, with e s 



In the case of 



Q q < Qf < u) m , the oscillating term can be neglected, 
leading to 

6 S = (5/2){uj m /n F ) 2 = 5Q m /(2n th + 1), (14) 

where we have defined the thermal occupation number 
nth through ksT / (huj m ) = nth + (1/2). For strong in- 
teraction Q. q 3> £If, the transcendental equation can be 
solved numerically, showing that 8 S > 1 is always valid. 
Maximally Entangled Mode. To gain insights into this en- 
tanglement, we apply the techniques in Ref. [28[ and de- 
compose outgoing fields into independent single modes by 
convoluting them with some weight functions namely 



6 l = (fS), [d i ,d}] = 26 ij , 



(15) 



which requires (fi\fj) 



If we define gn = 5R[/s] and 



g i2 = the single-mode quadratures will be 

X t ee (Oi + 6\)/V2 = J^dtga h - g i2 b 2 , (16) 
% ee (Oi - 6j)/(iV2) = f^dtga h + 9ll b 2 . (17) 



Different choices of weight function will generally give 
optical modes that have different strength of entangle- 
ment with the mechanical oscillator. The function of 
particular interest is the one that gives an effective opti- 
cal mode maximally entangled with the oscillator. Using 
the fact that logarithmic negativity is an entanglement 
monotone, the optimal weight function can be derived 
from the following constrained variational equation: 



{6Eff b /5g i )+mg i = Q (i = l,2), 



(18) 



where we have neglected unnecessary indices and /!& is 
Lagrange multiplier due to the constraint (f\f) = 1 
and Ej^ h quantifies entanglement in the subsystem con- 
sisting of the oscillator and the effective optical mode 
[i(0), p(Q), X, Y]. As it turns out, the optimal weight 
functions g^ 2 have the shape of decay oscillation with 
poles to given by the following polynomial equation 

[(« - uo m f + 7 2 J [(w + u m ? + lU + X = 0, (19) 

where parameter x is a functional of g± t2 and also depends 
on Q q and Qf- Therefore, the weight functions are 



Sfc(t) = A k e*» t coB(u g t + k ) (A = 1,2), 



(20) 



with 7 g and ui g being imaginary and real parts of tu. Ana- 
lytical solutions to parameters Ak ,u> gi ^ g and 9k requires 
exact expression of x in terms of gk, Q q and Q,p, which 
is rather complicated. Instead, we numerically optimize 
those parameters to maximize EJ^ h . 

Taking into account (/|/) = 1, A\ and A 2 can be re- 
duced to a single parameter Q, which is defined through 



At 



47 3 (7 g 2 +^)cos 2 [C + fc(^/2)] 
7 2 + ^ + 7 2 cos(26» fe ) + 7ff u g sin(2ft%) ' 



(21) 



From Eq. ([T 



7^. In addition, a local 



unitary transformation (rotation and squeezing) will not 
change the simplectic eigenvalue. Without loss of gener- 
ality, we can fix that B\ = ir/2 and 6 2 = 0. Therefore, 
only two parameters u g and £ need to be optimized. 

In the special case of weak-interaction and low- 
thermal-noisc limit (O g , Qp *C w TO ), the optimal £ op t is 
equal to 7r/4, which indicates A\ sa A 2 = 2^/j^ for a 
high-Q oscillator. Besides, as shown in the upper panel 



of Fig. El the optimal w° pt 



/(*) 



Hm e 



leading to 



(22) 



Therefore, the optimal weight function has the same 
shape as Stokes and Anti-Stokes sideband modes. This is 
similar to what been obtained in Refs. [l8, 28]; However, 
due to causality, the weight function here is defined in 
£ 2 [— oo,0] rather than C 2 [— oo,oo] which is essential for 
defining sideband modes. 

In the case of strong interaction and high thermal noise 
(Qq, Q-f > w m ), the optimal u> g deviates from u> m and de- 
pends on flp and f2 g , as shown in the lower panel of Fig. 
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FIG. 3: Logarithmic negativity E$ b as a function of quantity 
(uig — uj m ) / uj m in the weak-interaction and low-thermal-noise 
case (upper panel) and strong-interaction and high-thermal- 
noise case (lower panel). In the first case, we have chosen 
Q m = 10 3 , Q q /uj m — Q,F/tjJm = 2 x 10~ 2 . In the second case, 
Q m = 10 6 (independent of Q m for higher Qm), Q q /u m = 
50, n F /uj m = 20 and ( = tt/3. 



Conclusion. We have demonstrated that quantum entan- 
glement exists universally in system with a mechanical 
oscillator coupled to continuous optical fields. The en- 
tanglement measure — logarithmic negativity displays an 
elegant scaling which depends on the ratio between char- 
acteristic interaction and thermal-noise frequency. Such 
scaling should also apply in electromechanical systems 
whose dynamics are similar to what we have considered. 
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[3l More generally, the optimal Copt = tt/3 and u>° pt can 
be fitted by « (0.640,% + 0.57 fi 2 ,) 1 / 2 . Correspond- 
ingly, the logarithmic negativity can be approximated as 

E% h a (1/2) ln[l + (15. Sl 2 /(13. 0% + O 2 ))], (23) 

which again manifests universality of the entanglement. 
Therefore, as long as the optimal weight function is cho- 
sen, one can always recover quantum correlations be- 
tween the oscillator and the outgoing fields. 

In principle, by choosing a weight function orthogonal 
to the optimal one obtained above, one can derive next- 
order optimal mode. Repeating this procedure will gen- 
erate a complete spectrum of effective optical modes or- 
dered by Ej^ h , which is analogous to obtaining wavefunc- 
tions and corresponding energy levels with variational 
method. This not only helps to understand the full en- 
tanglement structure but also sheds light on experimen- 
tal verifications of such universal entanglement. Rather 
than trying to recover the infinite-dimension covariancc 
matrix in Eq. we can apply right weight functions 
to extract different effective optical modes and form low- 
dimension sub-systems. Take sub-system consisting of 
the oscillator and the maximally entangled optical mode 
for instance, 4x4 covariance matrix can be determined 
by measuring correlations among different quadratures. 
This can be achieved by using a local oscillator with time- 
dependent phase, which allows to probe both mechanical 
quadratures [2^] and those of the effective optical mode. 
For expample, a quadrature = X sm<^ + Ycos£ can 
be measured with the following local oscillator light: 

L(t) oc Li (t) cos ojot + L2 (t) sin uj^t (24) 

with L\(t) = <7i(t)cosC + <?2(t)sin£ and L 2 (t) = 
32 (i) cos £ — <7i(i)sinf. Synthesis of multiple measure- 
ments will recover the covariance matrix that we need to 
verify the entanglement. 
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